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ABSTRACT 

Aims. We present one-dimensional aperture synthesis imaging of the red supergiant Betelgeuse (a Ori) with VLTI/AMBER. We 
reconstructed for the first time one-dimensional images in the individual CO first overtone lines. Our aim is to probe the dynamics of 
the inhomogeneous atmosphere and its time variation. 

Methods. Betelgeuse was observed between 2.28 and 2.31 with VLTI/AMBER using the 16-32-48 m telescope configuration 
with a spectral resolution up to 12000 and an angular resolution of 9.8 mas. The good nearly one-dimensional uv coverage allows us 
to reconstruct one-dimensional projection images (i.e., one-dimensional projections of the object's two-dimensional intensity distri- 
butions). 

Results. The reconstructed one-dimensional projection images reveal that the star appears differently in the blue wing, line center, 
and red wing of the individual CO lines. The one-dimensional projection images in the blue wing and line center show a pronounced, 
asymmetrically extended component up to ~1.3 while those in the red wing do not show such a component. The observed one- 
dimensional projection images in the lines can be reasonably explained by a model in which the CO gas within a region more than 
half as large as the stellar size is moving slightly outward with 0-5 km s _I , while the gas in the remaining region is infalling fast 
with 20-30 km s . A comparison between the CO line AMBER data taken in 2008 and 2009 shows a significant time variation in the 
dynamics of the CO line-forming region in the photosphere and the outer atmosphere. In contrast to the line data, the reconstructed 
one-dimensional projection images in the continuum show only a slight deviation from a uniform disk or limb-darkened disk. We 
derive a uniform-disk diameter of 42.05 ± 0.05 mas and a power-law-type limb-darkened disk diameter of 42.49 ± 0.06 mas and a 
limb-darkening parameter of (9.7 ± 0.5) x 10~ 2 . This latter angular diameter leads to an effective temperature of 3690 ± 54 K for the 
continuum-forming layer. These diameters confirm that the near-IR size of Betelgeuse was nearly constant over the last 18 years, in 
marked contrast to the recently reported noticeable decrease in the mid-IR size. The continuum data taken in 2008 and 2009 reveal 
no or only marginal time variations, much smaller than the maximum variation predicted by the current three-dimensional convection 
simulations. 

Conclusions. Our two-epoch AMBER observations show that the outer atmosphere extending to ~ 1.3— 1 .4 R* is asymmetric and its 
dynamics is dominated by vigorous, inhomogeneous large-scale motions, whose overall nature changes drastically within one year. 
This is likely linked to the wind-driving mechanism in red supergiants. 

Key words, infrared: stars - techniques: interferometric - stars: supergiants - stars: late-type - stars: atmospheres - stars: individual: 
Betelgeuse 



1. Introduction 

Red supergiants (RSGs) experience slow, intensive mass loss up 
to 10~ 4 Mq yr _1 , which is very important for understanding the 
final fate of massive stars. For example, our poor understanding 
of the RSG mass loss makes it difficult to estimate the main- 
sequence mass range of the progenitors of the Type IIP super- 
novae, which are the most common type of core-collapse super- 
novae. The mass loss also plays a significant role in the chem- 
ical enrichment of galaxies. Despite this importance, there are 
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no satisfactory theories for the RSG mass loss at the moment, as 
stressed by Harper (20101. 

Studies of the outer atmosphere, where the winds are accel- 
erated, are a key to tackling this problem. The outer atmosphere 
of RSGs has complicated structures. The UV observations of 
the well-studied RSG Betelgeuse (a Ori, Ml-2Ia-Ibe) with the 
Hubble Space Telescope reveal that the hot (-6000-8000 K) 
chromospheric plasma is more than twice as extended as the 
photosphere (Gilliland & Dupree 19961. However, radio contin- 
uum observations with the Very Large Array show that much 
cooler (-4000-3000 K) gas extends to several stellar radii (Lim 
et al. 119981 1. suggesting that the hot chromospheric plasma and 
cooler gas coexist. IR spectroscopic and interferometric studies 
of a few bright RSGs also show the presence of dense H2O gas in 
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the outer atmosphere, the so-called "MOLsphere", extending to 
~ 1.3-2.0/?+ with column densities on the order of 10 20 cm -2 and 
tempera tures o f 1500-2000 K (e. g., Ts uji I2000al I2000bl [20061 
Ohnaka 120041 Perrin et al. 120041 l2D07l . The nature of the cool 
gas in the outer atmosphere of Betelgeuse has also recently been 
probed with mid-IR [Fe II] emission lines (Harper et al. 2009a). 
Possibly the chromospheric plasma with a small filling factor is 
embedded in more abundant, cooler gas. This inhomogeneous, 
multi-component nature of the outer atmosphere is considered 
to play a crucial role in driving mass outflows in RSGs. There is 
observational evidence for the asymmetric, inhomogeneous na- 
ture of the circumstellar material. For Betelgeuse, Kervella et al. 
(2009 ) found a very faint plume extending to ~6 R+ in the near- 
IR, while a millimeter CO map shows a blob at ~5" (~235 R*) 
away from the star (Harper et al. 2009b'. The emission of the 
CO fundamental lines near 4.6 fim shows that the circumstellar 
envelope is approximately spherical within 3" (140 R*) but with 
signatures of mildly clumpy structures (Smith et al. 2009 1. 

High-spectral and high-spatial resolution observations of 
strong IR molecular lines are ideal for probing the physical prop- 
erties of the outer atmosphere. The near-IR interferometric in- 
strument AMBER (Astronomical Multi-BEam combineR) at the 
Very Large Telescope Interferometer (VLTI) is well suited for 
this goal with its high spectral resolution up to 12000 and high 
spatial resolution down to 1-2 mas with the current maximum 
baseline of 130 m. In 2008, we observed Betelgeuse in the CO 
first overtone lines near 2.3 yum with AMBER (Ohnaka et al. 
2009 hereafter Paper I). The high spectral resolution of AMBER 
allowed us to detect salient signatures of inhomogeneities in the 
individual CO lines and to spatially resolve the gas motions in a 
stellar photosphere (and also MOLsphere) for the first time other 
than the Sun. 

However, in 2008, we obtained data only at six uv points, 
which are insufficient for image reconstruction. In order to ob- 
tain a more complete picture of the dynamics of the inhomoge- 
neous outer atmosphere, we carried out new AMBER observa- 
tions of Betelgeuse with a better uv coverage in 2009. In this 
second paper, we report on the first one-dimensional aperture 
synthesis imaging of Betelgeuse in the CO first overtone lines, 
as well as on time variation in the dynamics of the stellar at- 
mosphere in an interval of one year. The paper is structured as 
follows. The AMBER observations, data reduction, and image 
reconstruction are outlined in Sect. |2] We describe the results 
about the time variations as well as the one-dimensional image 
reconstruction in Sect. [3] The modeling of the velocity field pre- 
sented in Sect. [4] is followed by the discussion on the dynamics 
in the extended outer atmosphere (Sect. [5}- Conclusions are pre- 
sented in Sect. [6] 

2. Observations 

2.1. AMBER observations and data reduction 

AMBER (Petrov et al. l2007b is the near-IR (1.3—2.4 /^m) 
spectro-interferometric instrument at VLTI, which combines 
three 8.2 m Unit Telescopes (UTs) or 1.8 m Auxiliary 
Telescopes (ATs). AMBER measures the amplitude of the 
Fourier transform — the so-called visibility or visibility ampli- 
tude and two observables that contain information about the 
phase of the object's Fourier transform: differential phase (DP) 
and closure phase (CP). The DP roughly represents how the ob- 
ject's phase in a spectral feature deviates from that in the con- 
tinuum. Non-zero DP represents information about the photo- 
center shift in a spectral feature with respect to the continuum. 



The CP is the sum of the measured Fourier phases around a 
closed triangle of baselines (i.e., <p\2 + y?23 + f3i), not affected 
by the atmospheric turbulence. The CP is always zero or n for 
point-symmetric objects, and non-zero and non-7r CPs indicate 
an asymmetry in the object. Moreover the CP is important for 
aperture synthesis imaging in optical/IR interferometry. 

Betelgeuse was observed on 2009 January 5 and 6 with 
AMBER using three ATs in the E0-G0-H0 linear array config- 
uration with 16-32-48 m baselines (AMBER Guaranteed Time 
Observation, Program ID: 082.D-0280, PL: K. Ohnaka). As in 
Paper I, we used the A'-band high-resolution mode (HR_K) with 
a spectral resolution of 12000 covering wavelengths from 2.28 
to 2.31 //m to observe the strong CO first overtone lines near 
the (2,0) band head. Fringes could be detected on all three base- 
lines without the VLTI fringe tracker FINITO. We obtained a 
total of 54 data sets on two half nights. The data sets taken more 
than ~2 minutes apart were treated as separate data sets, be- 
cause Betelgeuse is strongly resolved and the visibility varies 
noticeably even for a slight change in the baseline length (see 
Fig. |2](. As shown in Fig. [1] the observed uv points align ap- 
proximately linearly at position angles of 73 ± 2°. This linear uv 
coverage allows us to sample the visibility function densely at 
this particular position angle and to reconstruct one-dimensional 
projection images as described in Sect. 12.21 Each data set con- 
sists of 500 frames (NDIT) with each frame taken with a detec- 
tor integration time (DIT) of 120 ms. Sirius (a CMa, A1V, K 
= -1.4) was observed for the calibration of the interferometric 
data of Betelgeuse. We adopted the same angular diameter of 
5.6 + 0.15 mas from Richichi & Percheron (2005) as adopted 
in Paper I. We only used the calibrator data sets obtained just 
before and after each data set on Betelgeuse. A summary of the 
observations is given in Table I A. fl 

We reduced our AMBER data with amdlib ver.2.20, which 
is based on the P2VM algorithm (Tatulli et al. 120071 ). Some data 
sets, particularly those measured on the longest baseline and/or 
near the CO band head, are too noisy for the analysis. Therefore, 
we improved the SNR by binning the entire raw data (object, 
dark, sky, and P2VM calibration data) in the spectral direction 
with a running box car function as described in Paper I. We used 
different binnings with these spectral resolutions: 

1. Spectral resolution = 12000 (no binning), 8000, and 4800 
for the data sets #16, #17, #18, and #45-#49 for comparison 
with the 2008 data (Sect.OandO- 

2. Spectral resolution = 6000 for all data sets for the image re- 
construction in the continuum and in the individual, isolated 
CO lines (Sect. [Hand El Fig. El). 

3. Spectral resolution = 1600 for all data sets for the image 
reconstruction near the CO band head at 2.294 /mi (Sect. 13.61 
Fig.©. 

We checked for a systematic difference in the calibrated vis- 
ibilities and differential/closure phases by taking the best 20%, 
50%, and 80% of all frames in terms of the fringe SNR (Tatulli 
et al. 120071) . The difference between the results obtained with the 
best 20% and 80% frames is typically -10%. We took the best 
20% of all frames for our final visibilities to avoid the system- 
atic bias due to the rapid atmosphere and to keep decent SNRs 
in the final results. On the other hand, the calibrated differen- 
tial/closure phases do not show this systematic dependence on 
the frame selection criterion. Therefore, we included the best 
80% of all frames for the final DPs and CPs. The errors of the 
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resulting visibilities, DPs, and CPs were estimated in the same 
manner as in Paper I. 

We removed telluric lines from the observed spectra of 
Betelgeuse as best as possible by using Sirius as a spectroscopic 
standard star. The telluric lines identified in the spectrum of 
Sirius were also used for wavelength calibration. As a template 
of the telluric lines, we convolved the atmospheric transmis- 
sion spectra measured at the Kitt Peak National Observatory^ 
to match the spectral resolutions of the data. The uncertainty 
in wavelength calibration is 2.0 x 10~ 5 jum (2.6 kms -1 ). We 
note that the uncertainties in wavelength calibration for the 2008 
data and 2006 data in Paper I were mistakenly overestimated. 
The correct uncertainties in the 2008 data and 2006 data are 
2.3 x 10~ 5 fim (3.0 km s" 1 ) and 1.8 x 10~ 4 (25.1 km s" 1 ), 
respectively. 

2.2. One-dimensional image reconstruction 

The good linear uv coverage along the position angle of 73° 
shown in Fig. Q] provides an opportunity to reconstruct the so- 
called one-dimensional projection image, which is obtained by 
integrating the object's two-dimensional intensity distribution 
along the direction perpendicular to the linear uv coverage on 
the sky (central slice theorem or Fourier slice theorem). In other 
words, this one-dimensional projection image represents the 
two-dimensional intensity distribution compressed or squashed 
onto the linear uv coverage on the sky. For example, the one- 
dimensional projection image of a uniform disk is a semi-circle 
(see also the two-dimensional image of a limb-darkened disk and 
its one-dimensional projection image shown in Figs. IB. lb and 
IB. lb ). The reconstruction of one-dimensional projection images 
was first proposed for radio interferometry by Bracewell (119 56b . 
Whereas the information in the direction perpendicular to the 
baseline vector is lost in one-dimensional projection images, 
they still provide model-independent information about the geo- 
metrical extent and asymmetry of the object. The reconstruction 
of one-dimensional projection images from IR interferometric 
data or lunar occultation data has been carried out (e.g., Navarro 
et al. [19901 Leinert et al. [19911 Tatebe et al. 120061 Chandler et 
al. 120071) . 

We used the MiRA package ver. O.9.S0 (Thiebaut et al. r2008l 
to reconstruct one-dimensional projection images at each spec- 
tral channel (details of our image reconstruction procedure are 
described in Appendix 151. We first carried out the image recon- 
struction using computer-simulated data to examine effects of 
the uv coverage and reconstruction parameters such as the ini- 
tial model, prior, and regularization scheme on the reconstructed 
images. These tests with simulated data are crucial for examin- 
ing the credibility of aperture synthesis imaging particularly for 
objects with complex structures. 

With appropriate reconstruction parameters determined from 
these tests, we attempted to reconstruct one-dimensional projec- 
tion images from the observed 162 visibility amplitudes and 54 
CPs. While this worked well for the continuum, the reconstruc- 
tion in the CO lines turned out to be very sensitive to the recon- 
struction parameters. For example, depending on the size of the 
uniform disk used as the initial model, the reconstructed one- 
dimensional projection image in the CO lines shows a faint re- 
gion on the eastern or western side. Therefore, we used the self- 
calibration technique, which has recently been successfully ap- 

2 http://www.eso.org/sci/facilities/paranal/instruments/isaac/tools/ 
spectra/atmos_S_K.fits 

3 http://www-obs.univ-lyonl.fr/labo/perso/eric.thiebaut/mira.html 



plied to AMBER data for the first time by Millour et al. (I2011I I. 
We added modifications to their technique to deal with some is- 
sue specific to our data of Betelgeuse as described in Appendix 
ICl This technique allows us to restore the phase of the complex 
Fourier transform of the object's intensity distribution from the 
DP measurements. Image reconstruction with the complex visi- 
bility (i.e., visibility amplitude and phase) removes the ambigu- 
ity of the solution derived with the visibility amplitude and CP 
alone. 



3. Results 

3.1. Continuum data 

To compare with the /f-band continuum visibilities from the 
2008 data that were derived from the binned data with a spectral 
resolution 4800, we also derived the visibilities from the 2009 
data binned with the same spectral resolution. As in Paper I, 
we selected continuum points shortward of the CO band head 
at 2.294 jjm. For each data set, we averaged the visibilities over 
the selected continuum points. We took the simple mean of the 
errors as the errors in the average continuum visibilities without 
reducing by -yWcont, where N cont is the number of the selected 
continuum points. The reason is that the measurement errors are 
dominated by the systematic error in the absolute visibility cali- 
bration and do not become smaller by the averaging. We applied 
this averaging to the 2008 data as well. Since the different con- 
tinuum spectral channels correspond to slightly different spatial 
frequencies, we also averaged the spatial frequencies from the 
selected continuum points. 

Figure [2] shows the A'-band continuum visibilities measured 
in 2008 and 2009 as a function of spatial frequency. The figure 
reveals that the nearly linear uv coverage shown in Fig.Q~]enabled 
us to sample the visibility function quite densely from the first 
to the fifth visibility lobe. Uniform-disk fitting to the 2009 data 
results in a diameter of 42.05 + 0.05 mas with a reduced^ 2 of 
3.8. Fitting with a power-law-type limb-darkend disk (Hestroffer 
et al. |1997| l results in a limb-darkened disk diameter of 42.49 ± 
0.06 mas and a limb-darkening parameter of (9.7 ± 0.5) x 10~ 2 
with a better reduced^ 2 of 2.5. While the reduced^ 2 value is 
still higher than 1, Fig.|2]shows that the deviation from the limb- 
darkened disk is not strong, as found for the 2008 data. Only at 
the highest spatial frequency (i.e., the smallest spatial scale) is 
the deviation noticeable, but the errors are also large there. 

The limb-darkened disk diameter derived from the /f-band 
continuum data and a bolometric flux of (111.67 ± 6.49) x 
10~ 13 W crrT 2 (Perrin et al. 2004) lead to an effective tempera- 
ture of 3690+54 K. We propose this value as an effective temper- 
ature of the continuum-forming layer, approximately free from 
the effects of molecular lines. Perrin et al. (20041 modeled K- 
broadband interferometric measurements of Betelgeuse with a 
continuum-forming blackbody sphere and an extended molecu- 
lar shell. They derived 3690 + 50 K for the continuum-forming 
sphere. This value excellently agrees with our effective tempera- 
ture of the continuum-forming layer. Our effective temperature is 
slightly higher than the 3600 ± 66 K recently derived by Haubois 
et al. (120091) from the the //-band observations with the Infrared 
Optical Telescope Array (IOTA), but both agree within the un- 
certainties. 

Figure [2] reveals that the continuum visibilities show no 
or only marginal time variations between 2008 (green dots) 
and 2009 (red and blue dots) within the measurement errors. 
We compare this observational result with the current three- 
dimensional convection simulation for RSGs by Chiavassa et 
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al. (2009). The visibility predicted for 2.2 fim (Fig. 18 of 
Chiavassa et al. 120091 1. which approximately samples the con- 
tinuum, shows maximum time variations of, for example, ±40% 
in the third lobe. However, the visibilities observed in the third 
lobe (at spatial frequencies of 66-68 arcsec -1 ) show no time 
variation within the error bars (3-5%), and these error bars 
are 8-13 times smaller than the predicted maximum variation. 
While it is not very likely that we observed Betelgeuse at two 
epochs when it accidentally showed the same, weak deviations 
from the limb-darkened disk, this possibility cannot be entirely 
excluded. However, it is also possible that the current three- 
dimensional convection simulation predicts too pronounced sur- 
face structures and time variations owing to the gray approxima- 
tion adopted for the radiative transfer, as Chiavassa et al. (2009) 
mention. The most direct test for three-dimensional convection 
simulations is to measure the amplitude of the temporal fluctua- 
tions in the visibility and closure phase as well as the time scale 
of fluctuations by long-term monitoring observations and com- 
pare these with the model predictions. AMBER observations at 
more epochs would be necessary to draw a definitive conclu- 
sion about whether or not Betelgeuse seen in the continuum 
indeed shows much weaker inhomogeneities and much smaller 
time variations than predicted by the current three-dimensional 
convection simulation. 

The deviation from the limb-darkened disk in the /f-band 
continuum visibilities is lower than that observed in the H broad- 
band by Haubois et al. (2009). Their measurements show devi- 
ations of the visibilities from the limb-darkened disk as high as 
80-120% already in the fourth lobe (converted from the squared 
visibilities plotted in their Fig. 4), where our /f-band continuum 
data still follow the limb-darkened disk within the measurement 
errors of 5-10% except for the data points near the visibility null 
at ~100 arcsec -1 . The cause of this difference is not yet clear, 
because of a number of differences between their observations 
and ours (e.g., differences in the observed wavelengths, spec- 
tral resolution, position angle coverage), //-band observations 
with higher spectral resolution and/or ZT-band observations with 
a wider position angle coverage are necessary to clarify this is- 
sue. 

3.2. Long-term behavior of the near-IR and mid-IR angular 
diameters 

Figure [3] shows the Zf-band uniform-disk diameters of 
Betelgeuse from the literature and the archival data summerized 
in Paper I, together with the 1 1 /mi uniform-disk diameter pre- 
sented in Townes et al. (2009 ), who found a noticeable decrease 
in the 11 fj.m size in the last 15 years, and the one-epoch mea- 
surement by Perrin et al. (2007 1. In marked contrast to the notice- 
able decrease in the 1 1 yum size, the ZT-band diameter has been 
quite stable for the last 18 years with only a possible, slight long- 
term decrease. These results can be qualitatively explained as 
follows. While the size of the star itself has been stable over the 
last 18 years, the temperature and densities and/or shape of the 
outer atmosphere have changed significantly (e.g., decrease in 
temperature and/or density). Because the mid-infrared apparent 
size is largely affected by the MOLsphere and by dust (Ohnaka 
120041 Verhoelst et al. 120061 Perrin et al. 120071 the changes in 
the outer atmosphere lead to a noticeable change in the 1 1 fim 
size. This interpretation has also been recently reached by Ravi 
et al. (120 101 ) based on the estimation of the surface temperature 
seen at 1 1 /um. On the other hand, the angular size measured 
with the Zf-broadband filter is only slightly affected, because the 
strong molecular bands of CO and H2O are present only at the 



short and long wavelength edge of the K band. Detailed model- 
ing of the mid-IR interferometric and spectroscopic data, which 
is necessary to quantitatively derive the change in the physical 
properties of the MOLsphere and dust, is beyond the scope of 
this paper and will be pursued in a forthcoming paper. 

3.3. Significant time variation in the CO first overtone lines 

Figure |4] shows a comparison between the data taken in 2008 
and 2009 for four representative CO lines. The results for the 
2009 data were obtained from the merged data of the data sets 
#16, #17, and #18, which were taken at uv points very close to 
one of the data sets obtained in 2008 (data set #1 in Paper I). 
The data sets #45—#49 were also taken at uv points very close 
to the same 2008 data set, and the results from these data sets 
agree well with those shown in Fig. [4] We used the binning with 
the same spectral resolution as applied to the data from 2008: 
spectral resolution of 12000 (no binning) for the 16 m baseline 
data, 8000 for the 32 m baseline data, and 4800 for the 48 m 
baseline data and CP, respectively. 

In marked contrast to the continuum data, Fig. |4]reveals sig- 
nificant time variations in the CO line visibilities. The visibility 
within each CO line obtained on the 16 m baseline in 2008 was 
characterized by the maxima in the blue wing and minima in the 
red wing (black line in Fig. |4j)). In the 2009 data, the visibil- 
ity on the 16m baseline does not show the maxima in the blue 
wing anymore (red line in Fig. |4j?) and is characterized only by 
the minima in the red wing. The visibilities on the 32 m and 
48 m baselines also show time variations, although the data on 
the 48 m baseline are noisy. 

Time variations are even clearer in the DPs and CPs. Non- 
zero DPs were not detected in the CO lines on the 16m baseline 
in 2008, but the 2009 data show clear non-zero DPs in the CO 
lines. On the other hand, the non-zero DPs on the 32m baseline 
obtained in 2009 are much weaker than those in the 2008 data. 
The DPs on the 48m baseline as well as the CPs measured in 
2009 also show significant time variations. The non-zero DPs 
and non-zero/non-7r CPs indicate the asymmetry of the CO-line- 
forming region in 2009, as found in 2008. 

The observed spectra also reveal changes in the line profiles. 
The lines observed in 2009 are redshifted by ~6 km s _1 com- 
pared to those observed in 2008, as shown in Fig.|4^. The spectra 
taken on 2009 Jan 5 and 6 agree very well, although they were 
calibrated independently. This confirms that the redshift of the 
CO lines in the 2009 data is real. All these results suggest that 
the dynamics in the atmosphere of Betelgeuse has changed in an 
interval of one year. 

3.4. One-dimensional projection images in the continuum 

The reconstructed one-dimensional projection image in the con- 
tinuum at 2.30662 fim is shown by the black line in Fig.|5](com- 
parison between the observed interferometric data and those 
from the reconstructed image is shown in Fig. ID.lt . Also plotted 
is the one-dimensional projection image of the limb-darkened 
disk with the angular diameter of 42.49 mas and the limb- 
darkening parameter of 0.097 (gray line, overlapping with the 
green line at angular distances between +15 and -15 mas). 
These one-dimensional projection images are already convolved 
with the Gaussian beam with a FWHM of 9.8 mas as described 
in Appendix [B] Figure [5] shows that the stellar surface is well 
resolved with the beam size of 1/4 of the diameter of the stel- 
lar disk. The one-dimensional projection images obtained at all 
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continuum spectral channels agree well within the uncertainty of 
the reconstruction (~1%). 

The overall deviation from the limb-darkened disk is small, 
~5% on the eastern side (position angle = 73°). Because inho- 
mogeneities in the direction perpendicular to the baseline vec- 
tor on the sky are smeared out in the one-dimensional projec- 
tion image, we estimated the upper limit on the strength of in- 
homogeneities using a uniform disk with one Gaussian-shaped 
dark spot. The //-band image of Betelgeuse shows spots with 
FWHMs up to 10-1 1 mas (Haubois et al. |2009l ), which we adopt 
for our model. For a given amplitude of the spot, we generated 
10000 models with random positions of the spot and counted 
the number of models whose one-dimensional projection im- 
age shows deviation from the uniform disk smaller than 5%. 
For spot amplitudes lower than 20% of the stellar disk inten- 
sity, more than 70% of the models show deviations compatible 
with the observations. However, the fraction of these models is 
13%, 3%, and <0.01% for spot amplitudes of 30%, 40%, and 
50% of the stellar disk intensity, respectively. Therefore, we es- 
timate the amplitude of the spot to be smaller than 20-30% of 
the stellar disk intensity with flux contributions of 1.6-2.4%. 

3.5. One-dimensional projection images in the individual CO 
lines 

Figure [5] shows the one-dimensional projection images recon- 
structed in the blue wing, line center, and red wing within the CO 
line centered at 2.3061 pm because of the two transitions (2,0) 
R(26) and R(75) with a spectral resolution of 6000 (compari- 
son between the observed interferometric quantities and those 
from the reconstructed images is shown in Fig. ID.lt . The one- 
dimensional projection images are normalized with the peak in- 
tensities but are not artificially registered with one another, be- 
cause the relative astrometry is preserved thanks to the restored 
visibility phase. The round or blunt shape of the images primar- 
ily results from the projection of the two-dimensional images 
onto the baseline vector on the sky. The level of the image re- 
construction noise is estimated to be ~1.5% from the strengths 
of the artifacts in the entire field of view used for the reconstruc- 
tion (256 mas). The one-dimensional projection images in other, 
isolated CO lines above ~2.3 pm agree well with those shown 
in Fig. [5] which adds fidelity to the image reconstruction. These 
one -dimensional projection images represent the imaging of the 
photosphere and MOLsphere of an RSG,for the first time, in the 
individual CO first overtone lines. 

Clearly, the one-dimensional projection images in the blue 
wing and line center are more extended than those in the con- 
tinuum, with the extension of 15% and 30% on the eastern and 
western side, respectively, when measured at the noise level of 
the image reconstruction (~1.5% of the peak intensities). The 
actual geometrical extension with intensities lower than the re- 
construction noise level could be even larger. These values are 
roughly consistent with the size of the MOLsphere probed with 
th e H 2 Q bands (e.g., Tsuii l2000all2006l Ohnaka 2004] Perrin et 
al. 2004, 2007). On the other hand, the red wing one-dimensional 
projection image shows only a slight deviation from the contin- 
uum without a trace of the extended component. Furthermore, 
the different appearance of the extended component in the blue 
and red wing suggests that the vigorous, inhomogeneous gas mo- 
tions are present not only in the photosphere extending to ~ 1 . 1 
R+ (see Sect. 4 in Paper I) but also in the layers extending to 
~1.3 This is because if the inhomogeneous gas motions were 
present only in the photosphere, the extended component would 
appear the same in the blue and red wing. 



3.6. One-dimensional projection images in the CO band 
head 

Figure [6] shows the one-dimensional projection images recon- 
structed at four different wavelengths in the CO (2,0) band head 
with a spectral resolution of 1600 (see Fig. lD.2l for a comparison 
between the observed interferometric quantities and those from 
the reconstructed images). The peak of the one-dimensional pro- 
jection image at the band head is shifted to the east with respect 
to that in the continuum and has geometrical extensions of 18% 
and 13% on the eastern and western side, respectively, when 
measured at the noise level of the image reconstruction. 

In none of the reconstructed images in the continuum, CO 
lines, and CO band head did we detect a feature corresponding 
to the faint plume reported by Kervella et al. (|2009), although 
their (single-dish) observations were carried out coincidentally 
almost simultaneously with our AMBER measurements. This is 
presumably because of the very faint nature of the plume. Its 
intensity is below 1% of the center of the stellar disk even at 
HI 3 pm, where it appears the most prominent, and the plume is 
even less pronounced at 2. 12 and 2.17 pm. Such a faint structure 
is below the noise of the image reconstruction from the present 
data (1.5% and 3% in Figs.[5]and[6] respectively). However, fu- 
ture observations with better accuracies could reveal the pres- 
ence of the plume in the CO lines, which would be useful for 
understanding the nature of the plume. Furthermore, the recon- 
structed images in the CO band head with a spectral resolution 
1600 can be compared with future medium-spectral resolution 
AMBER observations. 

4. Modeling of the velocity field 

We used our stellar patch model presented in Paper I to char- 
acterize the velocity field in the photosphere and MOLsphere. 
This model consists of two CO layers that represent the photo- 
sphere and MOLsphere. An inhomogeneous velocity field is rep- 
resented by a patch (or clump) of CO gas moving at some veloc- 
ity different from the CO gas in the remaining region. Because 
the reconstructed one-dimensional projection images do not al- 
low us to know the actual number and shape of the patches, we 
assumed only one patch in our modeling to keep the number of 
free parameters as small as possible. Furthermore, for the tem- 
perature, CO column density, and the radius of the two layers, we 
used the same parameters as derived in Paper I for the following 
reason. As discussed in Sect. [3] the CO line profiles observed 
in 2008 and 2009 show little time variation in the line depth and 
width except for the redshift. This implies that the physical prop- 
erties of the photosphere and MOLsphere, such as the density 
and temperature, may not have changed significantly, although 
there must have been temporal and spatial fluctuations. The pa- 
rameters adopted from Paper I are as follows: the inner CO layer 
is assumed to be located at 1 .05 R* with a temperature of 2250 K 
and a CO column density of 5 x 10 22 cm 2 , while the outer CO 
layer is assumed to be located at 1 .45 R* with a temperature of 
1800 K and a CO column density of 1 x 10 20 cnr 2 . We also 
adopted a microturbulent velocity of 5 km s for both layers 
as in Paper I. This means that we attempt to explain the one- 
dimensional projection images observed in 2009 by changes in 
the velocity field, as well as in the position and size of the patch. 
The wavelength scale of the model spectra was converted to the 
heliocentric frame assuming a heliocentric velocity of 20 km s _1 
(Huggins 1987; Huggins et al. [19941 . 

Figure Q shows the two-dimensional images, one- 
dimensional projection images, and the spectrum predicted 
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by the best-fit model for the same CO line as shown in Fig. 
The model is characterized by a large, off-centered, circular 
patch of CO gas, which dominates the upper half of the stellar 
disk (Figs. |7}3 and|7]i). The CO gas within this patch is moving 
outward with a velocity of 5 kms _I , while the gas in the 
remaining region is moving inward faster with 25 kms _I . 
Figures [7^ and|7f show that the line profile and the wavelength 
dependence of the observed one-dimensional projection images 
within the CO line are reasonably reproduced, although the 
difference between the images in the blue wing and red wing is 
somewhat too pronounced, and the line profile is weaker than 
the observed data. We found out that the models with a patch of 
CO gas moving slowly outward at 0-5 km s _1 with the gas in 
the remaining region downdrafting much faster at 20-30 km s _1 
can reproduce the observed one-dimensional projection images 
reasonably. 

The above model can also explain why the MOLsphere in 
the blue wing is much more pronounced than in the red wing. 
Firstly, the velocity field with the weak upwelling and strong 
downdrafting components causes the line center to be redshifted 
with respect to the stellar rest frame. This can be seen in Fig. [7k 
where the positions of the two transitions responsible for the ob- 
served line profile, (2,0) R(26) and R{75), are marked in the stel- 
lar rest frame. Since the contribution of the R(26) transition to 
the line profile is much larger than that of R(75) because of the 
much lower excitation potential of the former transition, the line 
profile is redshifted with respect to the position of R(26) in the 
stellar rest frame. In other words, the wavelength of the stellar 
rest frame is located in the blue wing of the line. Secondly, the 
strong extended CO emission is seen in the line of sight tangen- 
tial to the outer CO layer, because the column density along this 
line of sight is the largest. The radial velocity of the CO layer 
along such a line of sight is nearly zero, which means that the 
strong extended CO emission appears at the stellar rest frame. 
Because the stellar rest frame is located in the blue wing as ex- 
plained above, the extended CO emission is strong in the blue 
wing. The emission becomes nearly absent in the red wing, be- 
cause the velocity difference between the blue and red wing is 
much greater than the line width. 

The velocity field in 2009 is in contrast with that in 2008, 
which was characterized by the gas moving both outward and 
inward with velocities of 10-15 km s . Therefore, our AMBER 
observations at two epochs reveal a drastic change in the velocity 
field in the photosphere and MOLsphere within one year. 

5. Discussion 

The drastic change in the velocity field between 2008 and 2009 
sets an upper limit of one year on the time scale of the change of 
the MOLsphere. This allows us to estimate the upper limit of the 
radial spatial scale where the inhomogeneous gas motions are 
present. We assume that the upwelling patch (or clump) with 1 0- 
15 km s 1 detected at 1 .45 R* (radius of the MOLsphere) in 2008 
decelerated linearly with time over one year and corresponds to 
the patch slowly moving outward with 0-5 km s in 2009. Then 
the maximum radial distance reached by the gas patch is 0.24- 
0.47 R+. This means that the upwelling gas patch at 1.45 in 
2008 can reach 1.7-1.9 in 2009. Likewise, if we assume that 
the fast downdrafting patch with 20-30 km s detected in 2009 
accelerated inward linearly with time starting from km s _1 , it 
must have traveled 0.71-1.07 R+ in 1 year. This suggests that the 
fast downdrafting gas patch could have been located as far as at 
2.2-2.5 R+ in 2008 and could have fallen to 1.45 R+ in one year. 
Therefore, the vigorous gas motions can be present up to ~2 



These inhomogeneous gas motions in the extended atmo- 
sphere of Betelgeuse have also been detected by other obser- 
vations. Recently Harper et al. (2009a) have studied the dynam- 
ics of the cool extended outer atmosphere of Betelgeuse based 
on high-spectral resolution mid-IR observations of the [Fe II] 
lines at 17.94 and 24.53 yum. These [Fe II] lines form at ~1.6 
R+ (converted with the angular diameter of Betelgeuse derived 
here) in the cool extended outer atmosphere (see Fig. 8 of Harper 
et al. |2009al > with estimated excitation temperatures of 1520- 
1950 K. Therefore, the [Fe II] lines originate in the region sim- 
ilar to the MOLsphere where the CO first overtone lines form. 
The profiles of the [Fe II] lines indicate turbulent gas motions 
without signatures of significant outflows of ^;10 km s _I . This is 
consistent with the velocity fields derived from our two-epoch 
AMBER observations. Harper et al. (2009a) detected no signif- 
icant changes in the [Fe II] line profiles at three epochs over 14 
months. However, this may be because the changes in the veloc- 
ity field are smeared out in their spatially unresolved observa- 
tions. As can be seen in Fig. [4^, the CO line profiles observed 
with AMBER only show a low redshift, despite the remarkable 
change in the velocity field. 

Complex gas motions have been detected in the extended 
chromosphere of Betelgeuse as well. Lobel & Dupree (2001) 
present the modeling of the chromospheric velocity field up to 
~3 R+. Moreover, the velocity field changed from overall in- 
ward motions to outward motions within 0.5-1 year. Therefore, 
both the cool and hot components in the extended outer atmo- 
sphere are characterized by strongly temporally variable inho- 
mogeneous gas motions. 

The physical mechanism responsible for these vigorous mo- 
tions and their drastic change within one year is not yet clear, 
although it is likely related to the unknown wind-driving mech- 
anism. The convective energy flux is expected to be low in 
the MOLsphere, which extends to ~ 1.3-1. 4 This poses a 
problem for the interpretation of the detected gas motions in 
terms of convection. Other possible mechanisms include Alfven 
waves and pulsation. The recent detection of magnetic fields 
in Betelgeuse, albeit weak (~1 G), indicates that the prerequi- 
site for Alfven-driven-winds is available (Auriere et al. 1201 Oi l. 
Airapetian et al. (|2p00) show that Alfven waves can drive mass 
outflows from the chromosphere with the velocity and mass-loss 
rate in agreement with those observed for Betelgeuse. However, 
the effects of the Alfven-waves on the more dominant, cool outer 
atmosphere including the CO MOLsphere are not addressed. 
The MHD simulations of Suzuki ( 120071) for red giants show 
that the stellar winds are highly temporally variable and "struc- 
tured", in which hot (^ 10 4 K) gas bubbles are embedded in cool 
(~2000 K) gas (see, however, Airapetian et al. 120 101 for critical 
discussion). The radial velocity within ~10 also shows sig- 
nificant time variations from — 1-40 km s _1 (outward motions) 

to 40 km s _I (inward motions). This is compatible to the 

change in the velocity field detected by our AMBER observa- 
tions. However, the simulations of Suzuki (2007) were carried 
out for red giant stars, which are much less luminous (5 10 3 L0) 
compared to Betelgeuse (1.3 x 10 5 Lq, Harper et al. |2008 l. 
Extending the MHD simulations of Suzuki (2007) for more lu- 
minous stars, as well as the inclusion of the cool molecular com- 
ponent in the work of Airapetian et al. (2000 ), would be valuable 
for a comparison with the present and future AMBER observa- 
tions. 

Lobel (120101) infers that strong shock waves generated by 
convection in the photosphere that are propagating outward may 
carry the energy and momentum to accelerate the wind and heat 
the chromosphere. The qualitative similarity between the inho- 
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mogeneous velocity field in the chromosphere and in the pho- 
tosphere/MOLsphere may point toward this scenario. However, 
obviously, it is indispensable to map the dynamical structure of 
the cool outer atmosphere at various radii to clarify the driving 
mechanism of mass outflows in RSGs. 



6. Concluding remarks 

We have succeeded, for the first time, in one-dimensional aper- 
ture synthesis imaging of Betelgeuse in the individual CO first 
overtone lines, as well as in the continuum approximately free 
from molecular/atomic lines, with a spatial resolution of 9.8 mas 
and a spectral resolution of 6000 using VLTI/AMBER. 

The one-dimensional projection images in the CO lines re- 
constructed with the self-calibration technique, which restores 
the complex visibility using differential phase measurements, re- 
veal that the star appears different within the individual CO lines. 
The one-dimensional projection images in the blue wing and line 
center show a pronounced extended component up to 1.3 
while the images in the red wing follow that in the continuum 
without an extended component. Our image reconstruction rep- 
resents the first study to image the so-called MOLsphere of an 
RSG in the individual CO first overtone lines. Our modeling sug- 
gests that the dynamics in the photosphere and MOLsphere in 
2009 is characterized by strong downdrafts with 20-30 km s _I 
and slight outward motions with 0-5 km s _1 . This indicates a 
drastic change in the velocity field within one year from 2008, 
when the dynamics was characterized by both upwelling and 
downdrafting components with 10-15 km s _1 . 

On the other hand, the reconstructed one-dimensional pro- 
jection images in the A'-band continuum show only a small devi- 
ation of 5% from the limb-darkened disk with an angular diam- 
eter of 42.49 + 0.06 mas with a power-law-type limb-darkening 
parameter of (9.7 ± 0.5) x 10~ 2 . This limb-darkened disk diam- 
eter results in an effective temperature of 3690 + 54 K for the 
continuum-forming layer. The deviation from the limb-darkened 
disk in the one-dimensional projection images suggests that the 
amplitude of stellar spots is likely smaller than 20-30% of the 
intensity of the stellar disk. Furthermore, we detected no or only 
marginal time variation in the continuum visibility data within 
the measurement errors, much smaller than the maximum varia- 
tion predicted by the current three-dimensional convection sim- 
ulations. It cannot be entirely excluded that Betelgeuse showed 
unusually weak surface structures at the times of our AMBER 
observations just by chance. However, it is also possible that 
the current three-dimensional convection model for RSGs pre- 
dicts too strong surface structures in the continuum. A long-term 
monitoring to measure the amplitude of the time variations in the 
visibility and closure phase is indispensable for a definitive, sta- 
tistical test of three-dimensional convection simulations. 

The self-calibration imaging using differential phase has 
turned out to be very effective and necessary, despite the good 
linear uv coverage from the first to fifth visibility lobe. This sug- 
gests that the self-calibration technique may be even more nec- 
essary for two-dimensional imaging, where it is difficult to ob- 
tain a uv coverage as densely sampled as in our one-dimensional 
case. While the imaging of stellar surfaces is still challeng- 
ing (e.g., Creech-Eakman et al. 2010), our self-calibration one- 
dimensional imaging demonstrates a promising way to achieve 
that goal. 



reconstruction software MiRA publicly available. NSO/Kitt Peak FTS data on 
the Earth's telluric features were produced by NSF/NOAO. 
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Fig. 1. uv coverage of our AMBER observations of Betelgeuse. 
The dashed line represents the average position angle of 73°. 
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Fig. 2. Continuum visibilities of Betelgeuse averaged over the continuum spectral channels between 2.28 and 2.293 fim. The insets 
show enlarged views of the second, third, and fourth/fifth lobes. The solid and dashed lines represent the visibilities for a uniform 
disk with a diameter of 42.05 mas and for a limb-darkened disk with a diameter of 42.49 mas and a limb-darkening parameter of 
0.097 (power-law-type limb-darkened disk of Hestroffer 1997 ), respectively. The dotted lines represent the maximum range of the 
variations in the 2.22 jjm visibility due to time-dependent inhomogeneous surface structures predicted by the three-dimensional 
convection simulation of Chiavassa et al. ( 120091 1, who presents the model prediction up to 70 arcsec 1 . 
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Fig. 3. Uniform-disk diameters of Betelgeuse measured in the K 
band and at 1 1 fxm as a function of time. The black filled circles 
and the blue open circle represent the 11 fim diameters mea- 
sured by Townes et al. (2009) and Perrin et al. (2007i, respec- 
tively. The diamonds represent the Zf-band diameters from the 
foll owing references. IRMA: Dyck et al. (1 19921 1. IOTA: Perrin et 
al. d2004b . VINCI: Paper I. AMBER MR.K: AMBER medium- 
resolution (A/AA = 1500) data from Paper I. AMBER HR.K 1: 
AMBER high-resolution data from Paper I. AMBER HR_K 2: 
Present paper. The errors in the IOTA and AMBER measure- 
ments are smaller than the symbols. 
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Fig. 4. AMBER data in the CO first overtone lines obtained in 2008 and 2009. The black and red solid lines represent the data 
taken in 2008 and 2009, respectively. In the panels b-h, the scaled spectrum observed on 2009 Jan 5 is overplotted by the blue solid 
line to show the asymmetry of the visibilities and differential/closure phases within each line profile. The spectral resolutions given 
in the panels are the values of the binned data as described in Sect. 12.11 a: Normalized spectra. Two spectra derived from the data 
taken on 2009 Jan 5 and 6 are plotted by the red and green solid lines, respectively, while the 2008 spectrum is plotted by the black 
solid line, b-d: Visibilities observed on the E0-G0-16m, G0-H0-32m, and E0-H0-48m baselines, e: Closure phase, f-h: Differential 
phases observed on the E0-G0-16m, G0-H0-32m, and E0-H0-48m baselines. 
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Fig. 5. One-dimensional projection images re- 
constructed at four different wavelengths (blue 
wing, line center, red wing, and continuum) 
within the CO line due to the two transi- 
tions (2,0) R(26) and R(15) (spectral resolu- 
tion of 6000) are shown by the blue, green, 
red, and black solid lines, respectively. The 
one-dimensional projection image of the limb- 
darkened disk with the parameters derived in 
Sect. 13.11 is also shown by the gray solid line 
as a reference. The inset shows the observed 
line profile with the wavelengths of the im- 
ages marked by the filled circles with the corre- 
sponding colors. The one-dimensional projec- 
tion images are convolved with the Gaussian 
beam whose FWHM (9.8 mas) is shown by 
the thick solid line, and their absolute scale 
is normalized to the peak intensities. The ori- 
entation of the one-dimensional projection im- 
ages is also shown. The dashed line denotes the 
noise level of the image reconstruction. 
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Fig. 6. One-dimensional projection images re- 
constructed at four different wavelengths in 
the CO (2,0) band head (spectral resolution of 
1600) shown in the same manner as Fig. [5] 
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Fig. 7. Best-fit stellar patch model with an inhomogeneous velocity field, a: The filled diamonds represent the observed spectrum, 
while the black solid line represents the model spectrum. The filled circles mark the wavelengths of the two-dimensional images 
and one-dimensional projection images shown in the panels b-g. b-e: Two-dimensional model images. The orientation of the linear 
uv coverage is shown by the solid lines. The one-dimensional projection images are obtained by integrating the two-dimensional 
images in the direction shown by the dotted lines, which is perpendicular to the orientation of the linear uv coverage, f: The one- 
dimensional projection images predicted by the model at four wavelengths within the CO line, g: The observed one-dimensional 
projection images for the same CO line. 
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Appendix A: Summary of AMBER observations 

Our AMBER observations of Betelgeuse and the calibrator 
Sirius are summarized in Table IA.fl 

Appendix B: Image reconstruction of simulated 
data 

Aperture synthesis imaging from optical/IR interferometric data 
depends on a number of parameters used in the image recon- 
struction process, such as the initial model, as well as the reg- 
ularization scheme and prior, which represent a priori informa- 
tion about the object's intensity distribution. Therefore, it is im- 
portant to carry out the image reconstruction for simulated im- 
ages to derive the appropriate reconstruction parameters before 
we attempt the image reconstruction from observed data. From 
simulated images, we generate simulated interferometric data by 
sampling the visibilities and CPs at the same uv points as inter- 
ferometric observations. With the true images known for these 
simulated data, we can examine the appropriate reconstruction 
parameters that allow us to reconstruct the original images cor- 
rectly. 

Because MiRA is developed for two-dimensional image re- 
construction, we took the following approach for the reconstruc- 
tion of one-dimensional projection images: two-dimensional im- 
age reconstruction was carried out using an appropriate ini- 
tial model and regularization parameters as described below. 
The reconstructed two-dimensional image was convolved with 
the clean beam, which is represented by a two-dimensional 
Gaussian with a FWHM of /i/B max — 9.8 mas, where B max is the 
maximum baseline length of the data. The one-dimensional pro- 
jection image was obtained by integrating this convolved two- 
dimensional image in the direction perpendicular to the linear 
uv coverage. 

We generated two simulated images that represent possible 
surface patterns of Betelgeuse: a simple limb-darkened disk and 
a uniform disk with a bright spot, a dark spot, and an extended 
halo, as shown in Figs. IB. lb and IB. 2b , respectively. For both 
cases, the stellar angular diameter was set to be 42.5 mas, which 
is the limb-darkened disk diameter derived from all continuum 
visibilities measured in 2009. The visibilities and CPs were com- 
puted from the simulated images at the same uv points as our 
AMBER observations, using the program of one of the authors 
(K.-H. Hofmann). Noise was also added to the simulated visibil- 
ities and CPs to achieve SNRs similar to the AMBER data. We 
tested different initial models, priors, and regularization schemes 
to find out the appropriate parameter range to reconstruct the 
one-dimensional projection image of the simulated data cor- 
rectly. It turned out that uniform disks with angular diameters 
between 34 and 50 mas serve as good initial models. The prior 
used in the present work is a smoothed uniform disk described 
as 



where r is the radial coordinate in mas, and r v and e p define the 
size and the smoothness of the edge (e p — > corresponds to a 
uniform disk), respectively. The appropriate values for r p and e p 
were found to be 10 .. 15 (mas) and 2 .. 3 (mas), respectively. 
Therefore, we used six different parameter sets for the image re- 
construction of Betelgeuse by combining three diameters for the 
initial uniform-disk model (34, 42, and 50 mas) and two differ- 
ent parameter sets for the prior ((r p ,e p = (10, 2) and (15, 3)). 
The final images and their uncertainties were obtained by taking 



the average and standard deviation, respectively, from the results 
reconstructed with these six parameter sets. The regularization 
using the maximum entropy method turned out to be appropri- 
ate for our reconstruction. We started the reconstruction with a 
high degree of regularization (ju = 10 5 , see Thiebaut 120081 for 
the definition of ju) and reduced it gradually by a factor of 10 af- 
ter every 500 iterations until the reduced x 2 reaches ~ 1 or MiRA 
stops the iteration. These tests with the simulated data also con- 
firm the validity of our approach to reconstruct one-dimensional 
projection images using the MiRA software for two-dimensional 
image reconstruction. 

Appendix C: Self-calibration imaging with 
differential phase 

Because the principle of the self-calibration technique using DP 
measurements is described in detail in Millour et al. (201 Q, we 
mention the actual procedure only briefly. Then we describe the 
modification we added to this technique to deal with an issue 
specific to the AMBER data of Betelgeuse. 

The DP measured with AMBER at each uv point contains 
information on the phase of the complex visibility function and 
roughly represents the difference between the phase in a spec- 
tral feature and that in the continuum. However, two pieces of 
information are lost because of the atmospheric turbulence: the 
absolute phase offset and the linear phase gradient with respect 
to wavenumber. We can derive this lost phase offset and gradi- 
ent by a linear fit to the phase (as a function of wavenumber) 
from the reconstructed continuum images, if the image recon- 
struction in the continuum is reliable and not sensitive to the 
reconstruction parameters. This is indeed the case for our image 
reconstruction of Betelgeuse in the continuum, as discussed in 
Sect 13.41 Therefore, the phase in the CO lines can be restored 
from the continuum phase interpolated at the line spectral chan- 
nels and the DPs measured in the lines. The image reconstruction 
is carried out with the measured visibilities and CPs as well as 
the restored phase. This process can be iterated, but our exper- 
iments show that the reconstructed images do not change after 
the first iteration. 

We added the following modification to the technique pre- 
sented in Millour et al. (2011). When the phase offset and gra- 
dient are derived by a linear fit to the phase of the reconstructed 
images, we only use the continuum spectral channels below 

2.293 fim and those between the adjacent CO lines above 2.3 //m, 
instead of using the entire spectral channels, as Millour et al. 
( 1201 II) did. The reason for this selection of the spectral chan- 
nels is that the image reconstruction near the CO band head at 

2.294 fim is so uncertain owing to the poor SNR in the data 
binned with a spectral resolution of 6000 that the inclusion of 
the spectral channels near the band head in the linear fit hampers 
the reliable derivation of the phase offset and gradient. 

The inclusion of only the selected continuum channels has 
the following consequence. If we denote the continuum phase 
from the reconstructed continuum image at a given baseline and 
at the z'-th spectral channel as <p c (i), the phase at the z'-th spectral 
channel, tp(i), is restored as 

ip(i) = <p c (i) + DP(/), 

where DP(z) represents the differential phase at the z'-th spectral 
channel measured at the same baseline. At a continuum spectral 
channel denoted as z" c , the restored phase <p(i c ) should be equal to 
the phase from the reconstructed continuum image (p c (i c )- This is 
fulfilled if the measured DP in the continuum is zero. However, 
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Table A.l. Log of AMBER observations of Betelgeuse and the calibrator Sirius with the E0-G0-H0 (16-32-48m) baseline configu- 
ration. Seeing is in the visible. 
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(UTC) 


(m) 
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Betelgeuse 


2009 Jan 5 


2009 Jan 6 


1 


01:01:44 


22.24/11.12/33.36 


73 


0.76 


6.2 


500 


24 


00:47:33 


21.24/10.63/31.87 


73 


1.55 


3.8 


500 


2 


01:03:56 


22.45/11.23/33.68 


73 


0.77 


6.1 


500 


25 


00:49:46 


21.46/10.74/32.20 


73 


1.53 


4.0 


500 


3 


01:06:09 


22.66/11.33/33.99 


74 


0.77 


6.1 


500 


26 


00:51:59 


21.68/10.84/32.52 


73 


1.42 


4.7 


500 


4 


01:08:22 


22.86/11.44/34.30 


74 


0.84 


5.4 


500 


27 


00:54:12 


21.89/10.95/32.85 


73 


1.50 


4.7 


500 


5 


01:10:35 


23.07/11.54/34.61 


74 


0.89 


5.3 


500 


28 


00:56:25 


22.11/11.06/33.17 


73 


1.41 


5.0 


500 


6 


01:50:20 


26.41/13.21/39.62 


75 


0.58 


11.1 
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29 


00:58:37 
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73 
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4.8 
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7 


01:52:34 
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75 
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11.2 
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30 
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8.4 
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8 


01:54:46 
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75 


0.57 


11.6 
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31 


01:50:47 


26.73/13.37/40.11 


75 


1.10 


8.4 
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9 


01:56:59 


26.90/13.46/40.36 


75 


0.58 


12.2 


500 


32 


01:53:00 


26.89/13.45/40.35 


75 


1.01 


9.4 


500 


10 


01:59:12 


27.06/13.54/40.59 


75 


0.65 


11.2 
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33 


01:55:13 


27.05/13.53/40.59 
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9.5 
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11 


02:51:26 
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75 
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34 
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12 


02:53:39 
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75 
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75 
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16.4 
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13 


02:55:52 


30.31/15.16/45.47 


75 


0.91 


7.1 


500 


36 


02:35:07 


29.51/14.77/44.28 


75 


1.26 


14.2 


500 


14 


02:58:04 


30.40/15.21/45.61 


75 


1.02 


6.3 


500 


37 


02:37:20 


29.63/14.82/44.45 


75 


1.30 


11.8 


500 


15 


03:00:18 


30.49/15.25/45.75 


75 


0.98 


6.3 


500 


38 


02:39:33 


29.74/14.88/44.61 


75 


1.31 


11.8 


500 


16 


04:07:52 


31.98/16.00/47.98 


73 


1.19 


5.0 


1000 


39 


02:41:46 


29.84/14.93/44.77 


75 


1.45 


10.2 


500 


17 


04:12:08 


31.99/16.01/47.99 


73 


1.12 


5.3 


1000 


40 


03:20:21 


31.30/15.66/46.97 


74 


1.11 


17.8 


500 


18 


04:16:24 


31.99/16.00/47.99 


73 


1.09 


5.6 


1000 


41 


03:22:34 


31.36/15.69/47.06 


74 


1.02 


17.0 


500 


19 


05:00:17 


31.41/15.72/47.13 


71 


0.89 


12.4 


500 


42 


03:24:47 


31.42/15.72/47.14 


74 


0.99 


16.3 


500 


20 


05:02:29 


31.36/15.69/47.05 


71 


0.92 


12.4 


500 


43 


03:27:00 


31.48/15.75/47.22 


74 


1.01 


15.6 


500 


21 


05:04:42 


31.30/15.66/46.96 


71 


0.94 


12.1 


500 


44 


03:29:13 


31.53/15.77/47.30 


74 


0.90 


16.0 


500 


22 


05:06:55 


31.24/15.63/46.87 


71 


0.97 


13.0 


500 


45 


04:05:36 


31.98/16.00/47.99 


73 


0.97 




500 


23 


05:09:08 


31.17/15.60/46.77 


71 


0.92 


14.8 


500 


46 


04:07:49 


31.99/16.01/47.99 


73 


1.08 




500 
















47 


04:10:01 


31.99/16.01/48.00 


73 


1.16 
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48 


04:12:14 


31.99/16.01/47.99 


73 


1.11 




500 
















49 


04:14:27 


31.98/16.00/47.99 


73 


1.15 




500 
















50 


04:49:21 


31.57/15.80/47.37 


72 


1.40 




500 
















51 


04:51:34 


31.53/15.78/47.30 


71 


1.10 




500 
















52 


04:53:47 


31.48/15.75/47.23 


71 


1.02 




500 
















53 


04:56:00 


31.42/15.72/47.15 


71 


0.89 




500 
















54 


04:58:13 


31.37/15.70/47.06 


71 


0.92 




500 


Sirius 


2009 Jan. 5 


2009 Jan. 6 


CI 


01:28:18 


24.75/12.38/37.13 


50 


0.84 


5.9 


2500 


C5 


00:24:19 


20.99/10.50/31.50 


36 


1.32 


3.9 


2000 


C2 


02:17:35 


27.84/13.93/41.76 


59 


0.96 


8.5 


5000 


C6 


01:27:00 


24.92/12.47/37.38 


51 


1.34 


5.9 


2500 


C3 


03:23:41 


30.83/15.43/46.26 


66 


1.04 


6.0 


3000 


C7 


02:09:58 


27.57/13.79/41.36 


58 


1.02 


12.9 


2500 


C4 


04:36:43 


31.99/16.01/47.99 


73 


0.90 


7.5 


2500 


C8 


02:56:21 


29.91/14.97/44.88 


64 


2.01 


54.0 


2500 
















C9 


03:43:15 


31.46/15.74/47.20 


69 


0.81 


22.4 


2500 
















C10 


04:26:45 


31.99/16.01/47.99 


72 


1.56 


41.2 


2500 



the measured DPs show noticeable non-zero values in the con- 
tinuum, as exemplarily shown in Figs. IC. lb and lC.lb . The reason 
for the non-zero DPs in the continuum is that amdlib derives dif- 
ferential phase by a linear fit to the instantaneous phase at all 
spectral channels. Owing to the strong deviation of the phase in 
many CO lines from that in the continuum, this linear fit does 
not go through all the continuum points. Therefore, the non-zero 
DPs in the continuum spectral channels lead to a systematic error 
in the phase restored in the continuum, which affects the subse- 
quent image reconstruction. We found out that the continuum 
one-dimensional projection image reconstructed using the re- 
stored phase shows a systematic wavelength dependence from 
the shortest to the longest wavelength of the observed spectral 
range, which is not seen in the continuum images reconstructed 
from the visibilities and CPs alone. 

It is necessary to use the same spectral channels in the linear 
fit to the phase for the derivation of DP and for the derivation of 



the phase offset and gradient. Therefore, we refitted the DP from 
amdlib with a linear function (with respect to wavenumber) at 
the same continuum points as used for the derivation of the phase 
offset and gradient (dashed lines in Figs. IC.lfa and IC. lb ) and 
subtracted the fitted linear function from the observed DP. This 
procedure enforces the DP in the continuum spectral channels 
to zero within the measurement errors, as shown in Fig. IC.lt . 
The phase was restored using this "refitted" DP. The continuum 
one-dimensional projection images reconstructed using the refit- 
ted DPs do not show the aforementioned systematic wavelength 
dependence, which proves the validity of our procedure. 
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Fig. B.l. Image reconstruction of the simu- 
lated data for a limb-darkened disk with the 
parameters derived in Sect. 13.11 a: Original 
two-dimensional image of the simulated data. 
The solid line represents the orientation of 
the linear uv coverage, while the dotted line 
represents the orientation perpendicular to it. 
b: Two-dimensional image of the simulated 
data convolved with the Gaussian beam with 
a FWHM of 9.8 mas. c: Comparison between 
the original and reconstructed one-dimensional 
projection images before convolving with the 
Gaussian beam. The one-dimensional projec- 
tion images are obtained by integrating the 
two-dimensional images in the direction shown 
by the dotted lines in the panels a and b. 
d: Comparison between the original and re- 
constructed one-dimensional projection images 
convolved with the Gaussian beam with a 
FWHM of 9.8 mas. e: The filled circles and tri- 
angles represent the visibilities from the orig- 
inal simulated data and the reconstructed im- 
age, respectively, f: The filled circles and trian- 
gles represent the CPs from the original simu- 
lated data and the reconstructed image, respec- 
tively. The abscissa is the spatial frequency of 
the longest baseline of each data set. 
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Fig. B.2. Image reconstruction of the simulated 
data for a uniform disk with a bright spot, a 
dark spot, and a halo shown in the same manner 
as in Fig. IB. II 
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Appendix D: Fit to the interferometric data with the 
reconstructed images 



180 



Fi g tire s IDTI and IET21 show the fit to the observed interferometric 
data for the one-dimensional projection image reconstruction in 
the CO line and in the CO (2,0) band head shown in Figs. [5]and 
[6] respectively. 
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Fig. C.l. a: DP observed on the longest baseline in the data set 
#9 (Bp = 40.36 m) is plotted by the red solid line. The black line 
represents the scaled observed spectrum. The DP and spectrum 
are binned with a spectral resolution of 6000. The dashed line 
represents the linear fit at the selected continuum points as de- 
scribed in Appendix ICl DP = is shown by the dotted line, b: 
Enlarged view of the panel a for the CO lines. Note that the DP 
in the continuum points between the adjacent CO lines deviates 
from zero, c: DP after subtracting the linear fit to the continuum 
points as described in Appendix [C] is shown by the red solid 
line. The black line represents the scaled spectrum. The DP in 
the continuum points between the adjacent CO lines is now zero 
within the measurement errors. 
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Fig. D.l. Comparison between the observed interferometric data and those from the one-dimensional projection image reconstruc- 
tion for the CO line shown in Fig. [5] The first, second, third, and fourth columns show the comparison for the blue wing, line center, 
red wing, and continuum, respectively. The panels in the top row (a-d) show the observed CO line profile, and the filled circles 
denote the wavelength of the data shown in each column. In the remaining panels, the observed data are represented by the red 
circles, while the values from the image reconstruction are shown by the blue triangles. The reduced^- 2 values for the fit including 
the complex visibilities, squared visibilities, and CPs are given in the panels m-p. 
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Fig. D.2. Comparison between the observed interferometric data and those from the one-dimensional projection image recon- 
struction near the CO (2,0) band head shown in Fig. [6] The first, second, third, and fourth columns show the comparison for the 
continuum, blue side between the continuum and the band head, bottom of the band head, and red side of the band head, respectively. 
The panels in the top row (a-d) show the observed spectrum of the CO band head, and the filled circles denote the wavelength of the 
data shown in each column. In the remaining panels, the observed data are represented by the red circles, while the values from the 
image reconstruction are shown by the blue triangles. The reduced^- 2 values for the fit including the complex visibilities, squared 
visibilities, and CPs are given in the panels m-p. 



